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I. INTRODUCTION 

The theory of open quantum systems deals with the 
dynamics of quantum systems interacting with their sur- 
roundings. The most common approach to the descrip- 
tion of the time evolution of the system stems from the 
fact that the total system, i.e., subsystem plus environ- 
ment, is closed. Hence the density matrix, containing all 
the information on the state of the total system, satisfies 
the Liouville-Von Neumann Master Equation. However, 
due to the eventual complexity of the total system, this 
Master Equation usually cannot be solved neither ana- 
lytically nor numerically Moreover, the density matrix 
of the total system often contains much more information 
than what we actually need, since, by definition, we call 
'system' that part of the universe which is of interest for 
our study . In other words, we focus our attention only 
on the degrees of freedom of the system. Mathematically 
this corresponds to a trace over the environmental vari- 
ables leading to the reduced density matrix of the system 
ps- The dynamics of ps can be rather involved because 
of the effects of the interaction with the environment. In 
general, understanding these effects is not an easy task. 

In this paper we concentrate on a paradigmatic model 
of the theory of open systems, namely the damped har- 
monic oscillator. We consider a single harmonic oscillator 
interacting with a quantized environment modelled as an 
infinite chain of non-interacting oscillators. This model, 
also known as quantum Brownian motion (QBM) model, 
is central in many physical contexts, e.g., in quantum 
field theory , quantum optics 0, 0, an( A solid state 
physics |l| . The importance of the damped harmonic os- 
cillator model is also due to the fact that it is one of the 
few non-trivial systems for which an exact Master Equa- 



* Electronic address: sabrina@fisica.unipa.it 



tion for the reduced density matrix can be formulated 
and exactly solved ® Q, H S EH • For this rea- 
son it has been extensively studied both for fundamental 
and for applicative research. On the one hand, indeed, it 
allows to gain new insight in the process of decoherence, 
responsible for the appearance of a classical world from 
the quantum world |l3| . On the other hand, it is the 
key model for many types of new quantum technologies 
such as quantum computation EUj quan tum cryptogra- 
phy E3 1 Emd quantum teleportation |16| . 

In this paper we discuss a recent analytic approach 
[l2| aimed at solving the generalized Master Equation 
describing the reduced system dynamics in very general 
conditions. The method does not rely, indeed, neither 
on the Born-Markov approximation, leading to a coarse 
grained description of the dynamics, nor on the weak cou- 
pling assumption, valid for quasi-closed system. We will 
give an expression of the density matrix easily readable 
in physical terms since it is simply related to diffusion 
and dissipation coefficients. Moreover, we will use this 
analytic method to discuss some points usually mislead. 
In particular two different approximations both usually 
(and confusingly) referred to as Rotating Wave Approxi- 
mation (RWA) will be compared. We give examples clar- 
ifying the relationship between the following issues: non- 
Markovian dynamics, non-Lindblad type dynamics, and 
the positivity of the density matrix. 

The paper is organized as follows. In Sec. I we in- 
troduce the generalized Master Equation and its exact 
solution. In Sec. II we discuss the two types of RWAs 
and single out a class of observables not sensitive to the 
presence of rapidly oscillating terms. In Sec. Ill we study 
the short time non-Markovian dynamics of the mean en- 
ergy of the system and show the connection with Lind- 
blad or non-Lindblad type behavior. Finally, in Sec. IV 
conclusions are presented. 
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II. EXACT DYNAMICS 

Let us consider a harmonic oscillator of frequency ojq 
surrounded by a generic environment. The Hamiltonian 
H of the total system can be written as follows 



H = Hq + He + aXE, 



(1) 



where H = uj {P 2 + X 2 ) /2, H E and aXE are the sys- 
tem, environment and interaction Hamiltonians, respec- 
tively, and a is the dimensionless coupling constant. The 
interaction Hamiltonian considered here has a simple bi- 
linear form with position of the oscillator X and position 
environmental operator E. For the sake of simplicity 
we have written the previous expressions in terms of di- 
mensionless position and momentum operators and set 
H = 1. We denote the density matrix for the oscillator- 
environment system by p. 

Under the assumptions that: (i) at t = system and 
environment are uncorrelated, that is p(0) = /?s(0) £g> 
Pb(0), with pe the density matrix of the environment; (ii) 
the environment is stationary, that is [He, Pe(0)] = 0; 
(iii) the expectation value of the environmental operator 
E is equal to zero, that is Trg {Epe(0)} = (as for 
example in the case of a thermal reservoir ) , o ne can derive 
the following master equation for QBM jlq 



dps® 
dt 



= -Mo - 



A(i)(X 



S\2 



S-oS 



--r(t)(X 2 ) s + n(*)X S P E Jps(i). (2) 

We indicate with X s ( s ) and P s '( s ) the commutator 
(anticommutator) position and momentum operators re- 
spectively, and with Hq the commutator superoperator 
relative to the system Hamiltonian. Such a Master Equa- 
tion, obtained by using the time-convolutionless projec- 
tion operator technique 0, |2£j , is the superoperatorial 
version of the Hu-Paz-Zhang Master Equation Let 
us note, first of all, that the Master Equation is local 
in time, even if non-Mar kovian. This feature is typical of 
all the generalized Master Equations derived by using the 
time-convolutionless projection operator technique [2j,|21| 
or equivalent approaches such as the superoperatorial one 
presented in [12 ]. 

The time dependent coefficients appearing in the Mas- 
ter Equation are defined in terms of the noise and dissi- 
pation kernels [2| and contain all the information about 
the short time system-reservoir correlation. The coeffi- 
cient r(t) gives rise to a time dependent renormalization 
of the frequency of the oscillator. In the weak coupling 
limit one can show that this gives a negligible contribu- 
tion as far as the reservoir cut-off frequency remains finite 
0- The term proportional to 7(f) is a classical damp- 
ing term while the coefficients A(t) and H(t) are diffusive 
terms. 

The superoperatorial Master Equation can be ex- 
actly solved by using specific algebraic properties of the 
superoperators |T^ . The solution for the density matrix 



of the system is derived in terms of the Quantum Char- 
acteristic Function (QCF) Xt(x,p) at time t, defined by 



Ps 



W = 2^/ Xt(x,p) e-^ x - xP Uxdp. 



(3) 



It is worth noting that one of the advantages of the super- 
operatorial approach is the relative easiness in calculating 
the analytic expression for the mean values of observables 
of interest by using the relations: 



<n = (»r ( q-t^p) 



d" 



(4) 



x^p— 



Neglecting the frequency renormalization terms, the 
exact analytic expression for the time evolution of the 
QCF is [12 



Xt(x,p) = e- A ^ X o 



e- r ^x,e- r ^p 



(5) 



with xo QCF of the initial state of the system and 

x = cos(ojot)x + sm(ujot)p, 

p = — sm(ujot)x + cos(u>ot)p. (6) 

The quantity At(x,p) is a quadratic form in the position 
x and momentum p variables: 



t (x,p) = {x,p)A(t)( K X p ), 



(7) 



with 



A{t) 



e r{tl) R f (t, h)M(t)R(t, tx)dt x . (8) 



In this equation, the matrix R(t,ti) contains rapidly os- 
cillating terms and M(t) is given by 



A(t) -n(<)/2 \ 
-n(t)/2 o ) ' 



(9) 



with A(t) and H(t) diffusion coefficients present in the 
Master Equation @. Finally the time dependent term 
T(t) appearing in Eqs. JSJ and JSJ) is simply related to 
the dissipation coefficient j(t): 



T(t) 



7(ii)<fti 



(10) 



Eq. JSJ shows that the QCF is the product of an ex- 
ponential factor, depending on both the diffusion coef- 
ficients [A(t) and H(t)] and the dissipation coefficient 
[7(i)], and a transformed initial QCF. The exponential 
term accounts for energy dissipation and is independent 
of the initial state of the system. Information on the 
initial state is given by the second term of the product, 
the transformed initial QCF. Note that for long times 
Xo [e~ r (*'i, e~ r ^p] — > 1, and the system approaches, as 
one would expect, a thermal state at the reservoir tem- 
perature, whatever the initial state was. 



3 



III. ROTATING WAVE APPROXIMATIONS 

In this section we discuss the existence of two approx- 
imations often called with the very same name: RWA. 
Such a situation, of course, may cause misleading and, in 
some cases, can lead to an inaccurate description of the 
short time dynamics of a damped harmonic oscillator. 
A similar analysis on the RWA has been performed by 
Agarwal for the Master Equation describing spontaneous 
emission in a two-level system ^7) . 

Let us begin discussing what we will call the 'RWA 
performed after tracing over the environment'. Such an 
approximation consists in averaging over the rapidly os- 
cillating terms contained in the matrices R(t, t\), appear- 
ing in Eq. © (for the precise analytical calculation see 
Ref. ^3)- This approximation can be actually seen as a 
secular approximation, and therefore to avoid confusion 
we will call it hereafter with this name. Note that the mi- 
croscopic interaction Hamiltonian contains the so-called 
counter-rotating terms not conserving the unperturbed 
energy of the system and reservoir. The crucial point 
to stress is that the contribution of these terms is not 
washed out by the average over the rapidly oscillating 
terms, as we will show with a specific example later in 
this section. 

The second type of RWA is what we will call 'RWA 
performed before tracing over the environment' , or sim- 
ply the RWA. In this case the counter-rotating terms are 
already absent in the microscopic system-reservoir inter- 
action Hamiltonian which reads: 



H, 



E 



g n a 



*b n + h.c. 



(11) 



with b n annihilation operators of the reservoir harmonic 
oscillators, a = (X + iP) /V2 and a* = (X - iP) /V2. 
This interaction Hamiltonian is very common in Quan- 
tum Optics, but does not give an appropriate description 
of the dynamics for short times, since the virtual pro- 
cesses associated to the counter-rotating terms play an 
important role even for weak couplings. 

To illustrate better this point we consider the short 
time dynamics of a specific observable for the system: its 
mean energy Hq = ujQ(n + 1/2), with n quantum number 
operator. Ho belongs to a class of operators which do 
not depend on the rapidly oscillating terms averaged in 
the secular approximation ^3 ■ Hence for calculating the 
mean value of these operators we can use the simpler 
secular solution of the density matrix. The expression 
obtained in this way is exact. It has been shown [T^] that 
a sufficient condition to single out this class of operators 
is the following 



[(X s ) 2 -(P s ) 2 ]^ = 



X S P S A = 0, 



(12) 



with A generic observable of the system. 

Let us now compare the short time expressions of the 
mean quantum number n of the system obtained after 
performing the secular and RWA approximations. For 



simplicity we will consider the case of weak coupling be- 
tween system and a reservoir at T temperature. It is 
not difficult to prove that the short time non-Markovian 
expression of n(t), in the secular approximation, can be 
written as follows p"s| 



<n)(f ^a;- 1 ) 



2a l 



-4','UJ)! 2 ( n(w) + i ) (/.,' 



f 2 
2 ' 

where n(u) is the mean number of reservoir excitations 
at T temperature and g(uj) is the reservoir spectral den- 
sity with cut-off frequency uj c . For the considerations 
made earlier, this expression coincides with the exact 
mean value of n (see also ref. |22j). 

A similar calculation shows that, if we perform the 
'RWA before tracing over the environment', the short 
time expression of the mean number (n)awA{t -C w^ 1 ) is 
exactly half of (n)(t <C w" 1 ) as given by Eq. 1131) . Such a 
circumstance can be directly related to the fact that, in 
second order perturbation theory, the virtual processes 
due to the counter rotating terms combine to give rise 
to a real process. Thus, neglecting them, as it is done 
in Eq. Qll[l . amounts at neglecting one of the two chan- 
nels through which the system can exchange energy with 
the environment. For this reason the variation of energy 
predicted in the RWA is only half of the exact value. 



IV. NON-MARKOVIAN DYNAMICS 

In this section we will exploit the analytic solution pre- 
sented in Sec. II to clarify some common misbelief related 
to non-Markovian Master Equations and non-Markovian 
dynamics. In order to do this we consider again the spe- 
cific example presented in previous section, that is the 
dynamics of the mean energy of the system. In view of 
the considerations made before, we can simplify the dis- 
cussion concerning the dynamics of this observable by 
considering the solution of the Master Equation in the 
secular approximation 0, [2i| : 



dps A(t)+7(*) r t 

[aa'ps — 2a) psa + psaaJ](14) 



— = - - [a^aps - 2ap s a> + psa^a] 

A(t) - 7(*) 



As for Eq. we emphasize that this Master Equation, 
even if non-Markovian, does not contain reservoir mem- 
ory kernels usually associated to non-Markovian gener- 
alized Master Equations. In other words, Eq. (|14fl . as 
well as Eq. J5J, is local in time. By definition this means 
that such an equation of motion for the density matrix 
is characterized by the fact that the time derivative of 
Ps{t) only depends on the actual value of ps(t). On the 
contrary, in non-Markovian Master Equations involving 
memory kernels, the time derivative of ps{t) is related 
to values ps{t') of the density matrix at times t' < t. 
Dealing with a non-Markovian Master Equation local in 
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time is of course an advantage compared to a descrip- 
tion involving memory kernels, since the memory effects 
of the environment are incorporated in the time depen- 
dence of the coefficients. However, the locality of the 
Master Equations J5J and 114|) could seem hard to rec- 
oncile with the intuitive idea of the effects that a generic 
environment may produce. 

As noticed in Ref. |24|, it is the linearity of the mi- 
croscopic Hamiltonian model for QBM which forces the 
Master Equation to be local in time. Indeed, as a con- 
sequence of the linear microscopic system-reservoir in- 
teraction Hamiltonian, the density matrix of the reduced 
system is the solution of a linear integro-differential equa- 
tion (containing the memory kernel). The key point is 
that a function f(t) satisfying a linear integro-diffcrcntial 
equation, as the one satisfied by the density matrix, does 
not depend on its entire history. Its future behavior is 
uniquely determined by the Cauchy data, i.e., its initial 
value and the initial value of its derivative. For the case 
of a damped harmonic oscillator, the non-Markovian fea- 
tures are restricted by linearity to the 'memory'of the 
initial time instant [2J|. For this reason, the time de- 
pendent coefficient provide all the memory effects in the 
evolution of the density matrix. It is worth reminding 
that, as underlined by Paz and Zurek in |25| . 'perturba- 
tive Master Equations can always be shown to be local 
in time'. In the case of the damped harmonic oscillator 
this is true also for the exact Master Equation |25j . 

Another aspect worth stressing is related to the struc- 
ture of the Master Equations J2J) and its approximated 
form (|14fl . It is well known that, when the Markovian ap- 
proximation is performed, the Lindblad theorem ensures 
that the Master Equation describing the system dynam- 
ics can be put in the Lindblad form. In general, neither 
Eq. J2J) nor Eq. (|14|) . both non-Markovian, are (or can 
be recast) in Lindblad form. A question therefore may 
arise naturally: should we care about the positivity of the 
density matrix? The reason why we ask this question is 
related to Lindblad theorem, ensuring the positivity of a 
density matrix satisfying a Master Equation of Lindblad 
form. When the Master Equation is not in the Lindblad 
form, as it is in our case, there can be situations in which 
symptoms of unphysical approximations show up. For 
example, negative eigenvalues of the density matrix may 
appear. This is for instance the case of the Caldeira- 
Leggett model when some particular initial conditions, 
not consistent with the high temperature approximation, 
are assumed. 

In general, however, it is certainly not true that when- 
ever a Master Equation is not of Lindblad form then it 
does not preserve the positivity of the density matrix. 
A basic assumption of the theorem is indeed that the 
reduced system dynamics constitutes a dynamical semi- 
group 0- Thus, it can happen that the Master Equation 
is not in Lindblad form but conserves the positivity of 
the density matrix, while it violates the semigroup prop- 
erty. This is actually our case, as one can realize from the 
time dependence of the coefficients of the Master Equa- 



tions J3J and (|T4"|) . 

Now, let us have a closer look at the approximated 
Master Equation (|14|l . The form of this equation is sim- 
ilar to the Lindblad form, with the only difference that 
the coefficients appearing in the Master Equation 114|) 
are time dependent. For the sake of brevity we will de- 
note hereafter this Lindblad-type Master Equation with 
time dependent coefficients simply as Lindblad-type Mas- 
ter Equation whenever the time dependent coefficients 
are positive. Note, however, that Lindblad-type Master 
Equations, contrarily to Master Equations of Lindblad 
form (having constant coefficients), do not satisfy the 
semigroup property. The reason why we focus on the 
positivity of the time dependent coefficients of Eq. I|14|) 
is that the dynamics of the system depends crucially on 
their sign. One can define two different regions of the 
parameter space, the first one correspondent to the situ- 
ation in which the Master Equation coefficients are pos- 
itive (Lindblad type region) and the second one corre- 
spondent to the case in which the coefficients acquire, at 
some time instants, negative values (non-Lindblad type 
region). In Ref. 2&J a detailed study of the border 
between the Lindblad and non-Lindlad type regions is 
carried out and the reservoir parameters governing the 
passage from one region to the other one are singled out. 

The first and important property distinguishing the 
Lindblad from the non-Lindblad type regions is that in 
the first one, i.e. whenever the time dependent coeffi- 
cients of the Master Equation are positive, one may use 
the 'standard'Monte Carlo wave function method to un- 
ravel Eq. (|14|) |26j . and therefore to simulate the system 
dynamics. Such method is usually used to study numer- 
ically the time evolution of Markovian systems since, in 
order to apply it, one needs to cast the Master Equa- 
tion in the Lindblad form. Here we note that also when 
the Markovian approximation is not performed, if the 
Master Equation is of Lindblad type, the method can 
be applied. This is, for instance, the case of a damped 
harmonic oscillator interacting with a high temperature 
reservoir, as far as the reservoir cut-off frequency is bigger 
than the frequency of the system oscillator [23] . We note 
that up to recently, only quantum state diffusion unrav- 
elings were used for simulating the temporal behavior of 
a harmonic oscillator interacting with a non-Markovian 
quantized environment |27|. 

We now proceed to illustrate an example showing the 
short time dynamics of the system in the non-Lindblad 
type region, that is when the time dependent coefficients 
of the Master Equation are negative. To this aim we con- 
sider the case of the mean energy of the system oscillator 
when it interacts with an artificially engineered out of 
resonance reservoir. In order to study this example we 
first write down the density matrix solution in the secular 
approximation. 

This approximation washes out the contribution of the 
diffusive term H(t) in Eq. ©, which becomes 01 

A{t) = e~ r(t) / e r (' l) A(t 1 )^i = A r {t). (15) 
Jo 
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The simplified QCF thus reads 
Xt{x,p) 



e- r Wx,e- r Wp 



(16) 



By using Eq. Q one derives the following expression for 
the mean energy of the system ^1 

(H (t)) = e- T W(H Q (t = 0)> +LJ A r (t). (17) 

Let us now consider the case in which the system os- 
cillator is in its ground state, and the environment is in 
a thermal state at T temperature. We assume a Ohmic 
reservoir spectral density with Lorentz-Drude cut-off uj c . 
Moreover, we consider the case in which lo c < luo, that is 
the spectrum of the reservoir do not completely overlap 
with the frequency of the system oscillator. This is of 
course never the case in presence of a natural reservoir. 
However, it has been recently shown that for quasi-closed 
systems, as single trapped ions, it is possible to engineer 
different types of artificial reservoirs, and couple them in 
a controlled way to the ion motion [23 • 

Under the conditions stated above, for high reservoir 
temperature, the time dependent coefficient A(t) — j(t) 
shows an oscillatory behavior as depicted in Fig. ^ As 
one can see from this figure there exist intervals of time 
in correspondence of which A(t) — j(t) becomes nega- 
tive. Hence, in general, the Master Equation l|14f) is not 
of Lindblad-type. The density matrix, however, is al- 
ways positive and it is the semigroup property which is 
clearly violated, as one can see from Fig. |3 where the 
time evolution of the mean quantum number operator 
(n) (heating function) is plotted. The signature of the 
non-Lindblad type dynamics of the reduced system is 
given by the presence of oscillations in its mean energy. 
These oscillations are due to virtual exchanges of energy 
between the system and the reservoir. As we have shown 
in references [3 an d H3 > these virtual processes are ab- 
sent whenever the Master Equation (|14f> is of Lindblad- 
type. In Figure[21the analytic result given by Eq. I|17fl is 
compared with the simulation using the Non-Markovian 
Monte Carlo method 0, El ■ In this case indeed, as we 
have previously discussed, the 'standard 'Monte Carlo 
wave function method cannot be applied. 



FIG. 2: Short time behavior of the heating function for dimen- 
sionless coupling constant au>o = 0.1Hz, u>q = 10 7 Hz, r = 0.1, and 
T = 300K. We compare the analytic solution and the Monte Carlo 
simulation with 10 7 histories. 



V. CONCLUSIONS 

In this paper we have used the exact solution of a 
damped harmonic oscillator to discuss some relevant 
issues of open system dynamics often leading to confusion 
and/or misunderstandings. 

We have focussed our attention on the so called RWA 
and we have pointed out the existence of two different ap- 
proximations called with the very same name. In other 
words we have stressed the differences in what we call 
the l RWA performed before or after tracing over the en- 
vironment. The RWA performed before tracing over the 
environment consists in neglecting the counter-rotating 
terms in the microscopic Hamiltonian describing the cou- 
pling between system and environment. The RWA per- 
formed after tracing over the environment is more pre- 
cisely a secular approximation, consisting in an average 
over rapidly oscillating terms, but does not wash out the 
effect of the counter-rotating terms present in the cou- 
pling Hamiltonian. By considering a specific example we 
show how, in the short time non-Markovian regime, the 
virtual processes due to the counter-rotating terms (ne- 
glected in the RWA) give a significant contribution even 
for weak coupling, and thus need to be taken into ac- 
count. 

We have pointed out that non-Markovian Master 
Equations do not necessarily contain memory kernels. 
This concept, already claimed by many authors (see for 
example 0, |25j) is still - surprisingly - often seen in a 
rather sceptical way. We have also discussed a related 
point which often rises doubts: the conservation of the 
positivity of the density matrix when the Master Equa- 
tion ruling its dynamics is not in Lindblad form. 

In order to study the dynamical properties of non- 
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Markovian Master Equations, we have introduced two 
subclasses of Master Equations which are not in the Lind- 
blad form: the Lindblad and non-Lindblad type, and we 
have stressed their difference. We have presented a new 
example showing that there exist conditions under which 
the short time dynamics of a damped oscillator is gov- 
erned by a non-Lindblad type Master Equation whose 
solution is always positive. In this case the semigroup 
property for the reduced dynamics which is violated, and 
the non Lindblad-type dynamics shows up in the exis- 
tence of virtual energy exchanges between system and 
reservoir. 

We think that this paper contributes in clarifying some 
aspects, often mislead, of the non-Markovian dynamics 
of a damped harmonic. This is important also because 
very recently the potential interest in non-Markovian 



reservoirs for quantum information processing has been 
demonstrated (3fl | and a non-Markovian description of 
quantum computing, showing the limits of the Marko- 
vian approach, has been presented |3l) . 
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